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ABSTRACT 

We expound the impact of extra sterile species on the ultra high energy neutrino fluxes in neutrino 
telescopes. We use three types of well-known flux ratios and compare the values of these flux ratios 
in presence of sterile neutrinos, with those predicted by deviation from the tribimaximal mixing 
scheme. We show that in the upcoming neutrino telescopes, its easy to confuse between the 
signature of sterile neutrinos with that of the deviation from tribimaximal mixing. We also show 
that if the measured flux ratios acquire a value well outside the range predicted by the standard 
scenario with three active neutrinos only, it might be possible to tell the presence of extra sterile 
neutrinos by observing ultra high energy neutrinos in future neutrino telescopes. 
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1 Introduction 



Cumulative effort through a series of experiments spanning more than four decades, has helped 
augment our understanding of the properties of neutrinos. That neutrinos have mass, mixing and 
hence flavor oscillations, has been firmly established. The picture emerging from the combined 
results of solar [1], atmospheric [2], reactor [3, 4], and accelerator [5, 6] neutrino experiments, is 
seen to be consistent with Ara^ = 8 x 10~ 5 eV 2 and Am§ x = 2.5 x 10~ 3 eV 2 , 1 and the so-called 
"tribimaximal mixing" [7] pattern for the neutrinos [8, 9]. In the tribimaximal (TBM) mixing 
scheme the mixing matrix has the form 



In this scenario, the atmospheric neutrino mixing angle #23 is maximal (sin 2 6*23 = 0.5), the so-called 
Chooz mixing angle 13 = 0, while the solar neutrino mixing angle is such that sin 2 12 = 1/3, all 
consistent with the current global data set. Maximal #23 and zero #13 can be easily obtained by 
imposing discrete family symmetries such as the ji — r symmetry [10] or the L M — L T symmetry 
[11]. While n — t symmetry does not predict the neutrino masses and the solar mixing angle, 
a phenomenologically viable scenario which satisfies all neutrino data can be obtained within an 
approximate L M — L T scheme [11]. 

The TBM mixing scheme will be put to stringent test with the next generation experiments 
involving neutrinos from the sun, atmosphere, reactors and accelerators. The mixing angle sin 2 9 12 
can be measured to less than 16% accuracy at 3a using solar neutrinos and to 6% accuracy from 
a future ~ 60 km baseline SPMIN (Survival Probability MINimum) reactor experiment [12]. 
Another prospect for accurate determination of sin 2 12 is to dope the Super-Kamiokande or even 
megaton water detectors with gadolinium [13]. The small and hitherto undetermined mixing 
angle #13 will be probed in the up-coming reactor experiments [15], as well as in accelerator based 
experiments using beams from conventional sources [16], beams produced by decay of accelerated 
radioactive ions stored in rings ( "Beta-Beams" ) or beam from decay of accelerated muons stored 
in rings ("Neutrino Factory") [18]. In principle, one of the best determination of #13 could be 
possible from the neutrino signal of a future galactic supernova [19]. The deviation of #23 from 
maximality and the sign of -D 2 3 = 0.5 — sin 2 #23 can be experimentally checked in atmospheric 
neutrino experiments [14]. 

A very well known and interesting feature of the TBM mixing scheme arises in the context of 
the flavor ratios of ultra high energy neutrinos arriving on earth. Ultra high energy neutrinos are 
created mainly through decay of high energy pions produced in pp and J97 collisions. Hence their 
relative flavor content at the source is expected to be 






(2) 



1 We follow the convention where Am? = mf — 

1 j 1 
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where 0° (a = e, fi, r) are the fluxes at the source. Under the TBM scheme, due to the inherent 
ji — t symmetry, #23 = 45° and 9 13 = 0, and therefore one obtains the flux ratio at earth as [20] 



Next generation km 3 neutrino telescopes such as the IceCube in Antarctica [21], Km3NET in 
the Mediterranean [22] are being especially designed and built to detect ultra high energy neutrinos 
coming from astrophysical sources. A lot of interest has been recently generated on the potential 
of using the observed flavor ratios of the ultra high energy neutrinos in neutrino telescopes in 
deciphering the predictions or deviation from tribimaximal mixing [23]. The key idea is that the 
flavor ratios are predicted to be one, as given by Eq. (3), only if 9i 3 = and #23 = 45°, a prediction 
of TBM mixing. If the mixing matrix was to deviate slightly from the TBM prediction through 
a change in #23 or #13 , we would see a difference in the flavor ratios, which could be used to pin 
down the extent of this deviation. 

Here we study how the flavor ratios of ultra high energy neutrinos get affected if we have extra 
sterile neutrinos mixed with the three active ones. Sterile neutrinos have been the subject of 
much discussion recently, following the much-awaited MiniBooNE results [24]. The MiniBooNE 
experiment was designed to specifically test the observed signal at the LSND experiment [25], 
which can be explained most convincingly by neutrino oscillations. The MiniBooNE data seem 
to contradict the LSND signal, as they do not see the kind of electron excess predicted by the 
oscillation explanation of the LSND data sample. Since the claimed flavor oscillation observed by 
LSND demands a Am 2 ~ eV 2 , it cannot be accommodated along with the solar and atmospheric 
neutrino observations within a three-generation framework. Three separate Am 2 scales can be 
possible if we have at least 4 neutrino states. Since the number of light active neutrinos are 
restricted to three from the decay width of the Z bosons at LEP, the extra neutrinos have to 
be sterile. The most economical scenario, with only one extra sterile neutrino leads to the so- 
called 2+2 and 3+1 mass spectra for the neutrinos [26]. The 2+2 spectrum is heavily disfavored 
by the solar and atmospheric neutrino data [9]. The 3+1 mass spectrum though comparatively 
less disfavored before the MiniBooNE results, suffers from a tension in explaining simultaneously 
the observed oscillation signal in the LSND experiment and the null signals in the other short 
baseline experiments [27]. The full analysis of all short baseline data including the MiniBooNE 
results further disfavors the 3+1 scheme [28]. Adding two sterile neutrinos turns out to be the 
next option to reconcile the LSND observations with the rest of the world neutrino data. In 
this so-called 3+2 scheme [29] the tension between the LSND and MiniBooNE data is reduced if 
one allows for CP violation, and this mass spectra is found to be consistent with global neutrino 
oscillation results [28]. 

In this paper we look at the signatures of extra sterile neutrinos on the observed flavor ratios 
of ultra high energy neutrinos in neutrino telescopes. We look at the impact of sterile neutrinos 
on the flavor ratios vis-a-vis the impact due to deviation from TBM mixing. Firstly, for small 
active-sterile mixing, we elucidate the confusion between the two cases and identify the range of 
the mixing angles between the active and sterile neutrinos which would give flavor ratios close to 
the ones predicted by deviated TBM mixing scenarios. Next we show that large values of hitherto 
unconstrained mixing angles between the active and sterile neutrino species could lead to extreme 
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values for flavor ratios, providing smoking gun signal for the existence of light sterile neutrinos 
which are heavily mixed with the active neutrinos. Thus signal from neutrino telescopes could 
provide an independent check of the results from the short baseline experiments in general and 
LSND and MiniBooNE in particular. 

We present all results in this paper in the framework of 3+1 scenario just for simplicity. It 
is absolutely straightforward to extend our results and conclusions to the 3+2 scenario. Here we 
have chosen not to work with the 3+2 mass spectra because that would entail many extra mixing 
angles, making the results look complicated. Therefore for the sake of illustration, we allow only 
two active sterile mixing angles to be non-zero in this work. In section 2 we give a brief discussion 
of the ultra high energy neutrino fluxes and the corresponding flux/flavor ratios. In section 3 the 
prediction from TBM mixing and the impact of deviation from TBM mixing in the framework of 
just three active neutrinos is reviewed. Section 4 gives the oscillation probabilities for three active 
and one extra sterile neutrino, while section 5 gives the corresponding flux ratios. In section 6 we 
present our numerical results. We end in section 7 with our conclusions. 

2 Ultra high energy neutrino fluxes 

Ultra high energy neutrinos are predicted from a number of astrophysical sources. The fireball 
model which has been put forward as the most plausible source for gamma ray bursts, are expected 
to produce ultra high energy neutrinos as well. Protons moving in the jets of the fireball are 
accelerated to very high energies. These highly accelerated protons undergo pp (or pn) and 
collisions, producing pions. The pions produce neutrinos through their decay channel 

7r + — ► /i + + , followed by /i + — > e + + v e + (4) 

and 

7r~ — > /i~ + , followed by /i~ — ► + u e + . (5) 

It is expected that the pry reaction will predominantly produce 7r + , while both ir + and tt~ are 
expected from the pp and pn collisions. Therefore, for the pp and pn collisions we expect the 
flavor ratio 

{< : <Pl e : < : <t>% : < : < : } = {l : 1 : 2 : 2 : : o} , (6) 

where <^P Va and <pp a are the neutrino and antineutrino fluxes respectively, of flavor a. For pry 
collisions on the other hand, we expect 

{< : <Pl e : < : < : < : < : } = {l : : 1 : 1 : : o} . (7) 

Though apparently it looks that the flux ratio for the two dominant channels of ultra high en- 
ergy neutrinos are different, the neutrino telescopes will not have sensitivity to the charge of the 
resultant leptons and hence in general are not expected to be able to distinguish the neutrino of 
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a given flavor from its antineutrino. Only when the Glashow resonance (P e e — > W~ — > anything) 
is the detection channel, V e is unambiguously detected and it is possible to measure the v e and 
v e signal separately. However, this occurs only over a small energy window E = m^/2m e = 6.3 
PeV, where m w and m e are the W and e mass respectively. In what follows, we will consider the 
sum of the neutrino and antineutrino signal at the neutrino telescope and in general refer to the 
flavor ratio as 



= 1:2:0, (8) 



where (fi a = (p Ua + <p Pa . If the ultra high energy neutrinos come from a neutron source, then the 
flux ratio is expected to be 

< : 4>l e : < : < ■ <-t : < : } {() : 1 : : : : ()} . (9) 

while for a muon damped source where the secondary muon is trapped and hence does not decay, 
it would be 

{<:<:<:<:<:<:} = {0:0:l:l:0:0}. (10) 

However, in this paper we will assume that only pions are the sources of ultra high energy neutrinos 
and work with the flavor ratio given in Eq. (8). 

Since the absolute flux predictions for the ultra high energy neutrinos could be uncertain by 
a huge amount, it is better to work with ratios of the fluxes. In this paper we exemplify the 
predictions of deviation from TBM mixing, and the impact of an additional sterile neutrino, using 
the following flux ratios: 

Re = a ^ a ' R » = a t4 ' Rt = a "7~a ' ( n ) 

M + (f) r (p e + (p T (f) e + (f)^ 

where e , M and r respectively are the z/ e + z/ e , i^ + P^ and z/ r + z/ T fluxes at earth after oscillations. 
Written in terms of the oscillation probabilities: 

p _ -Pee + 2 P efl _ 2 Pjjjx + Pen „ _ 2 P^ T + P eT 

e OP _|_P _lP _i_ P ' M P +9P A- P A- IP ' T P A- IP A- IP A- P 

where we have used P efJi = P^ e , since there is no chance of CP violation here. For the case for 
three active neutrinos only the flux ratios are given by 

p \ ^ i^ejx ~ Per) p _ 2 — (P e /j + 2 Pfi T ) „ -Per + 2 P^ T ^ o \ 

e "2-(P e ,-P er ) ' l + (P eM + 2P, r ) ' ^-3-(P er + 2P, T ) ' Udj 

where P Q/ 3 is the z/ a — > i/g oscillation probability discussed in the following sections. For the case 
where we have one extra sterile neutrino in addition to the three active ones, the flux ratios can 
be written as 

j-y 1 "I - (-Pe^t Per) Pes 2 (-Pe/i 2P^ r ) 2P^ S Per "I - 2 P^r 

e O / D DA O D ' ^ 1 i / D i Q D \ D' r 



(Pe/i Per) 2P jtts 1 -f- (Pe^i "I - 2P /ttT ) Pes 3 (Per "I - 2P^ T ) Pes 2P jUS 

where P es and P MS are the z/ e — > z/ s and z/ M — > z/ s oscillation probabilities respectively. However, we 
prefer to express the flux ratios in the sterile case by the most general form given by Eq. (12). 
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3 Oscillation probabilities and flux ratios in TBM mixing 



For neutrinos coming from astrophysical sources the oscillation probabilities can be written in 
general as 

-Pa/3 = Sa/3- ^UaiUpiUajUpj 

= YH u -i\ 2 W, (14) 

i 

where i,j and a, (5 are the generation indices of the mass basis and flavor basis respectively. 
The oscillatory terms, sm 2 (Amf j L/AE), have averaged out to 1/2 and the hence the probabilities 
depend only on the elements of the mixing matrix. The neutrino mixing matrix for three active 
neutrinos can be parameterized as [31] 



/ C12C13 S12C13 Si^e lS \ 



U 



—S12C23 — Ci2>S'23S'i3e* C12C23 — <S'i2<S'23*S'i3e l S23C13 
V S12S23 — C12C23 Si3e lS —C12S23 — Si2C23Si3e lS C23C13 J 



(15) 



where = cos 9ij and Sij = sin fly , and 5 the CP phase. Under the tribimaximal mixing ansatz, 
^12 = V3j 5*23 = 1/2 an d Sf 3 = 0, hence Eq. (15) reduces to the form given by Eq. (1). The 
oscillation probabilities therefore have an extremely simple form and are given as 



5 pTBM ]_ pTBM _ 2 pTBM _ % pTBM _ J_ pTBM _ ]_ 

g' W* lg ' e M - Q , er -g, M r rr "jg- 



The effect of the fi — r symmetry in-built in TBM mixing is clearly reflected by P^ M = P^ M = 
pTBM anc ^ pT BM — P^ T BM in Eq. (16). Therefore, the flux ratios defined in the previous section 



are 



pTBM _ pTBM _ pTBM _ }_ ^yj^ 

Any deviation from TBM mixing in the three active neutrino case will be reflected in these flux 
ratios, which will deviate from 1/2. We calculate the probabilities using Eqs. (14) and (15) for 
the case where mixing angles are different from those predicted by TBM mixing and present the 
results for the flux ratios using Eq. (13). 



4 Oscillation probabilities with sterile neutrinos 

In what follows, we will work in a framework where we include one extra sterile neutrino in addition 
to the three active ones and parameterize the 4x4 mixing matrix as 

U a = R(634)R(624)R(023)R(0u)R(0l3)R(0l2) , (18) 

where R{9ij) are the rotation matrices and 8ij the mixing angle. Note that we have put all phases 
to zero in Eq. (18). In general for the 3+1 scenario there are 3 CP violating Dirac phases. 
However, in this paper we have neglected the effect of the CP violating phases for simplicity. 
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4.1 TBM mixing with sterile neutrinos 

If we assume that the mixing angles #12, $13 and #23 follow the same values as in TBM mixing, 
i.e., sin 2 # 12 = 1/3, sin 2 9 13 = and sin 2 #23 = 1/2, and we further assume that the mixing angle 
614 = 0, then the 4x4 mixing matrix U s is given by 



U 8 = 



_ V3 

|[C34 + S24S34] — \[\[Cz4 + S24S34] 

v \Zi"^ 34 ^ 24 ~ ^ 34 -l ~ \[\[Cz±S24 — 5*34] 

Using Eq. (14) we get the oscillation probabilities: 



2 rr2 
e2 u el ) 



1 - 2U^U t 
5 

9 ' 

-2C/ e2 [/ el C/ /l2 ^ 1 , 



-a 2 
9 



24 ' 



\ 

S24 

C24S34 

_ V^^ 34 ^ 24 ^ 34 -l ^24C34 J 



\[Czi — S24534] 



Pw = 1 



1 - 2C 2 + -C* 



18 



24 5 



2^3 {u% + t^) - 2f/ 2 4 (l - f/ 2 4 ) , 



(19) 



(20) 



(21) 



(22) 



p = 



P — 

1 /IT 



-2U e2 U el U T2 U 



Tl ) 



- (C34 + S24S34) 2 , 

y 



2C/ M2 C/ M lC/r2C/rl " 2[/ /l3 C/ r3 (^2^2 + U^U t1 ) + ^ 
— (^C^C^ + 25C 2 4 S , 24S'|4 — 4C2 4 C34S , 245 , 34) • 



18 



(23) 



(24) 



4.2 Deviation from TBM mixing and sterile neutrinos 

If the mixing angles #12, #13 and #23 are different from that predicted by the exact TBM ansatz, 
then one has to use the expression for the full U s given by Eq. (18). For simplicity we use just 
one case where the mixing in the 3x3 active sector could deviate from TBM, viz, the case where 
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the mixing angle #23 is non-maximal. For this case the matrix U s is given by 



U. 



( 

\[\[SlzCzi + 6*235*24 1S34 ] - 
V V3 [623634 5*24 — 523S34] - 

and the oscillation probabilities are 



623624 



H/S23634 - l 23^2 I J 



|f6236: 



23O34O24 



623 5*24 S; 
S24 — S23 5; 



23*J34 





5*23624 
[623634 — 5*235*245*34] 
-[634523524 + 623534] 



\ 

524 

624534 
624634 J 



,(25) 



5/9, 



p — 

1 e/i 



MM 



"23 w 24 ) 



P — 

J P.T 



P » T ~ 36^ 



^ r< 2 n 2 
9 

1 ^ „ 4 ^4 or 1 ^ c 2 ^ 2 

1 — -0 23 2 4 — Z(_/ 24 D23t_/. 

4 2 
- (634523 + 524534623) 

1 



23 



9 q2 

ZD24O24 , 



(26) 
(27) 

(28) 

(29) 



24 



566 2 2 36 3 2 4 5 2 2 3 + 5 24 {(57 - 862 23 + 764 2 3)5245 3 2 4 + 2(-2 + 7622 3 )522352 3 4}](30) 



where 62^ = cos 29^ and S2ij 



sin26ij. In this paper we will mainly discuss the deviation 



from TBM mixing by changing sin 2 #23 from its maximal value. The TBM mixing ansatz is also 
violated if either sin 2 #13 7^ or sin 2 9 12 7^ 1/3. We will present one plot where we show the flux 
ratios as a function of sin 2 We will see that the effect of deviation of the flux ratios from the 



TBM predicted value of 1/2 due to sin #13 is very small compared to the change due to sin #23 . 
We therefore keep sin 2 13 fixed at zero for the results predicted by the sterile case in this paper. 
2 Therefore for all results with sterile neutrinos, we will keep sin 2 13 = and sin 2 9± 2 = 1/3 for 
simplicity and allow only sin 2 #23 to deviate from maximality. 



5 Flux ratios with sterile neutrinos 

As discussed before, flux ratios for ultra high energy neutrinos are the best model independent 
probes for understanding neutrino properties. If we consider the case where the only deviation for 
TBM mixing comes from #23 being non-maximal, while sin 2 9 12 = 1/3 and #13 = 0, and continue 
to work under the approximation that #14 = 0, then the flux ratios R e , i? M and R T are given as 

_ 5 + 86 2 4 6 2 3 

e JJdenom ' ^ ' 

where, 

j^denom _ _|_ \c\^C\ z — 86336^ — 3662 4 5 2 3 623 — 365^6^ + 4(63 4 523 + 52453462s) 2 

— 8634623(634523 + 53452462s) 2 + 36624523623(623634 — 523524534) (524534623 + 523634) 
+ 365 2 4 5 3 2 4 6 2 2 4) . (32) 



2 Since sin 2 6*13 = the effect of the normal CP phase is also completely absent. 
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where, 

R denom _ g _|_ gC% 4 C% 3 + 4(C34<S23 + S24 5*34 C23) 2 + ^^^^ — 8Cf 4 Cf 3 (C34S23 + S24 5*34 C23) 2 

+ 18C , 2 4 S'223 (C34C23 — 5*23 5*34 )(S24>S34C23 + S23C34) , (34) 



and 



R num 

23 + I2C23C24 — 8C24C23 — 9C| 4 S'223 — 9 5224 



where, 



R num _ 4((7 34( S' 23 -\- S34S24C23) 2 + 9S'| 4 S'224 — 8^24^23 (^W^S + 5*34 5*24 C23) 2 

+ 18C , 2 4 S'223 (C34C23 — 5 , 23>S , 24>S , 34) (S24S34C23 + S23C34) . (36) 

5.1 Flux ratios when #24 = 

If we work in a further simplified scenario where the only non-zero sterile mixing angle is the angle 
#34 and put both 6*14 = and #24 = 0, then the flux ratios are given as 

Re = 5 + 8C 23 ( 37 ) 

22 — 4S , 2 3 S , 34 — 8C 23 — 3QS23C23 + 28S , 2 3 C23C34 

R = 18 + 4C 2 3 - 8C& - 9£2l 3 

5 + 8Cf 3 + 4C 3 45'23 + 71S223C34 



4C* 5*, + 28C2 2 3^3C3 2 4 ,, p v 

" " 23 -24^^ + 4^ ' 



5.2 Flux ratios when 6*34 = 

If instead of putting 6*24 = 0, we put #34 = along with 6*14 = then the flux ratios are 

Re = 5 + 8C| 4 C| 3 ^ 

22 — 4Cf 3 — 36C24S24 + 32C24CI3 — 28C24S24C23 — 36C24CI3 

R _ 18 + 4Cf 4 Cf 3 — 8C| 3 C| 4 — 9C| 4 5'22 3 — 95*2^4 
9 — 28CI4C23 — 4Cf 3 + 36CI4CI3 
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23 + I2CI3C24 — 8C24C23 — 9(7245*223 — 95224 

Since the constraints from the short baseline experiments restrict # 24 to very small values, we can 
expand the flux ratios in a Taylor series and keep only the first order terms in sin 2 #24, giving 



5 + 8C% 4 C% 3 
22 — 8CI3 — 365|4 + I2C235I4 + 105223524 



R ° I "^24^23 / 43 N 

1L e — 00 orl o o/?o2 , in/i4 o2 , 1 n or>2 o2 ' V °/ 



^ „ (18 - 32C& + 28C| 3 ) - 45 2 2 4 (9 - 17C| + 14C|) ^ 
9 — 28C| 4 C| 3 — 4C| 3 + 36CI4CI3 

i2 ~ 45| 3 + 28C| 4 C| 3 5| 3 

r (23 - 24C 2 3 + 28C 2 4 3 ) - 45 2 2 4 (9 - 15C 2 3 + 14C 2 4 3 ) ' 1 ' 

6 Results 

The mixing angles for both active-active and active-sterile mixing are constrained by the results 
of the neutrino oscillation experiments [1, 2, 3, 4, 5, 6, 25, 27]. The current 3a limits on 12 , 6> 2 3 
and 9 13 are [8, 9] 

0.25 < sin 2 9 12 < 0.39 , (46) 
sin 2 2fl 23 > 0.9 , (47) 

sin 2 9 13 < 0.044 . (48) 

For the parameterization given by Eq. (18) which we use for the mixing matrix corresponding 
to the 3+1 scenario, the mixing angles governing the active-sterile mixing are 9 U , 9 2 a and 6* 34 . 
The mixing angles 14 and #24 are severely constrained by LSND [25] and other short baseline 
data [24, 27]. We refer the reader to [28] for the most up-to-date analysis of the short baseline 
experiments. In what follows, we will present all our results for 

sin 2 9 14 = (49) 

and for very small values of sin 2 #24 . Both these mixing angles are severely constrained by the 
current data. However, the mixing angle #34 does not appear in either P ee or P efl which are probed 
by most of the short baseline experiments like Bugey, CDHS, KARMEN, LSND and MiniBooNE. 
This mixing angle appears in the P^ T and P er channels. These channels were probed by the Chorus 
and NOMAD experiments at CERN. However, in the mass square difference range relevant for 
LSND, the constraint on the mixing angle 6> 34 are extremely weak. We will therefore allow sin 2 #34 
to take all possible values. We reiterate that the 3+1 mass spectrum for the neutrinos is now 
comprehensively disfavored after the release of the MiniBooNE results. Our choice of still using 
it as an exemplary case stems purely from keeping the discussion simple. Our choice of taking 
9 14 = is also due to the same reason. In fact, we would have a similar situation even if we had 
taken the 3+2 mass scheme and kept all angles except #24 and #34 non-zero. 
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Figure 1: The flux ratios R e , and R T versus sin 2 #23, for three active neutrinos. The yellow 
(light) bands corresponds to allowed ranges of the flux ratios corresponding to all currently 3a 
allowed values of sin 2 612 and sin 2 6*13 and where cos 5 is kept fixed as 1. The deep green (dark) 
bands show the corresponding allowed values where sin 2 #12 and sin 2 #13 vary in their current 3a 
allowed ranges and cos<5 is allowed to vary between —1 and 1. The upper (lower) boundaries of 
the yellow and green bands almost coincide with each other in the left (middle and right) panels. 



6.1 Three active neutrinos and deviation from TBM mixing 

In Fig. 1 we show the flux ratios as a function of sin 2 #23 for three active neutrinos. The blue 
dashed lines show the case where sin 2 #12 and sin 2 #13 are kept fixed at their TBM values of 
1/3 and respectively. We see that for this case, R e varies between [0.437-0.569] as sin 2 #23 
changes from [0.4-0.6]. Correspondingly, R^ and R T vary between [0.488-0.558] and [0.446-0.511] 
respectively. We note that R e depends sharply on sin 2 #23 and decreases linearly as the value of 
sin 2 023 increases. Both R^ and R T show a non-linear increase as sin 2 #23 increases. R^ increases 
faster when sin 2 623 > 0.5, while R T shows a sharper dependence when sin 2 #23 < 0.5. This means 
that R e is a good probe of deviation of 6*23 from maximality and its octant for all values of sin 2 #23, 
while i? M is good only for sin 2 6*23 > 0.5 and R T for sin 2 #23 < 0.5. If we allow even sin 2 #12 and 
sin 2 #13 to deviate from their TBM mixing values and vary freely within their current 3cr ranges 
keeping 6 = fixed, then we obtain the range of flux ratios shown by the yellow (light) bands in 
Fig. 1. If the CP phase S is also allowed to take all possible values then the results obtained are 
shown by the green (dark) bands in Fig. 1. We see that the allowed range of R e , R^ and R T for 
the currently allowed 3a range of sin 2 823 between [0.34-0.66] are respectively 

0.378 < R e < 0.645 , 
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Figure 2: The flux ratios R e , and R T versus sin 2 #13, for three active neutrinos and when 
sin 2 9 12 = 1/3, sin 2 ^23 = 1/2. The black dotted lines corresponds to cos 5 = 1, the red solid lines 
to cos 5 = and the green dashed lines to cos 5 = — 1 . 



0.476 < Rf, < 0.644 , 

0.379 < R T < 0.524 . (50) 

where we allow sin 2 9%2 and sin 2 #13 to vary within their current 3a allowed values and 5 over its 
full range. If we had three active neutrinos only, since we know that sin 2 2#23 > 0.9 at 3cr, the 
flux ratios must be restricted within the bounds given by Eq. (50). If we measure a flux ratio 
outside these ranges in the future neutrino telescopes, then that would be a certain indication of 
new physics and in particular as we will see in the next subsection, of presence of extra sterile 
neutrinos mixed with the active ones. 

In Fig. 2 we show the flux ratios as a function of sin 2 #13, keeping sin 2 9 12 and sin 2 #23 at their 
TBM values. The black dotted lines in this figure show the case for cos S = 1, the red solid lines 
for cos 5 = and green dashed line for cos 5 = —1. A comparison of Fig. 2 with Fig. 1 shows that 
the impact of changing sin 2 #23 from its TBM value of 1/2 is much larger on the flux ratios than 
changing sin 2 # 13 from its TBM value of 0. Therefore, in the rest of the paper we only show the 
effect of deviation from TBM mixing by changing sin 2 #23 from 1/2 and keep sin 2 #13 fixed at 0. 

6.2 Impact of sterile neutrinos 

For the sake of illustration we begin by presenting results for the case where the active-sterile 
mixing angles sin 2 9u = and sin 2 #24 = 0. For the active-active mixing angles we take sin 2 9 12 
and sin 2 ^13 at their TBM values of 1/3 and respectively and allow only sin 2 #23 to vary on 
both sides of its TBM value of 1/2. In other words, deviation from TBM mixing is caused by 
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Figure 3: Contour bands of the flux ratios R e , and R T in the sin 2 #23 — sin 2 #34 plane. The 
color scheme corresponds to flux ratios increasing with red, orange, yellow, green and so on. The 
R e bands increase from 0.44 to 1.1 in steps of 0.066, the R^ bands increase from 0.49 to 1.18 in 
steps of 0.069, while the R e bands increase from 0.00 to 0.51 in steps of 0.051. We have kept 
#14 = = 824 and all other mixing angles at their TBM mixing values. 




Figure 4: Panel (a) shows the contours in the sin 2 2#23 — sin 2 #34 plane which satisfy the corre- 
sponding condition given by Eq. (51). Each point on a given line gives the value of sin 2 # 34 which 
along with TBM in active part would give the same flux ratio as three active neutrinos but with 
that value for sin 2 2#23- Panel (b) gives the corresponding contours in the sin 2 2623 — sin 2 624 plane, 
showing the condition Eq. (55). 
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a non- maximal #23, while presence of sterile neutrino and its impact comes through the non-zero 
mixing angle #34. 

In Fig. 3 we show the contour bands of the flux ratios R e , and R T in the sin 2 # 2 3-sin 2 9 34 
plane. The color scheme corresponds to flux ratios increasing with red, orange, yellow, green 
and so on. The R e bands increase from 0.44 (lower right-hand corner) to 1.1 (upper left-hand 
corner) in steps of 0.066, the R^ bands increase from 0.49 (lower left-hand corner) to 1.18 (upper 
right-hand corner) in steps of 0.069, while the R T bands increase from 0.00 (top) to 0.51 (bottom) 
in steps of 0.051. As noted in the previous subsection, we see from the figure that R e has a 
very sharp dependence on sin 2 # 23 and its value decreases as sin 2 823 increases. i? M shows a milder 
dependence and its value increases with sin 2 623, while R T is almost independent of the value of 
sin 2 623 for sin 2 623 > 0.5 and shows an extremely mild dependence for sin 2 623 < 0.5. All flux ratios 
depend strongly on the value of sin 2 634 . The flux ratios R e and R^ increase as sin 834 increases, 
while R T decreases. This can be seen also from Eqs. (37), (38) and (39), where sin 2 634 decreases 
the denominators of R e and R^, while it decreases the numerator of R T . The figure also shows 
that the dependence of R T on sin 2 634 is linear, while that of R e and R^ are non-linear. This is 
corroborated by the Eqs. (39), (37) and (38) respectively. The most important thing we can see 
from this figure is the following: the same value of R e , R^ and R T can be obtained for different set 
of values of sin 2 623 and sin 2 6 34 . This gives us degenerate solutions in the sin 2 6 l 23 -sin 2 6> 34 plane. 
This will therefore lead to confusion 3 between the case where there are only three active neutrinos 
with deviation from TBM mixing and the case where TBM holds in the three active regime but 
there are extra sterile neutrinos which are mixed with the active neutrinos. This degeneracy is 
obviously more pronounced when the correlation between sin 2 623 and sin 2 634 is greater. Therefore, 
the largest effect is expected in R e and smallest in R T . 

To further illustrate the confusion between deviation from TBM mixing and presence of sterile 
neutrinos we present Fig. 4(a), which shows in the sin 2 2# 23 -sin 2 634 plane the contours for the 
condition 

^(sin 2 # 23 ^ 0.5, sin 2 9 M = 0) = ^(sin 2 9 23 = 0.5, sin 2 9 U ^ 0) , (51) 

where a = e, ji or r. The L.H.S of Eq. (51) gives the flux ratios predicted by deviation from 
TBM mixing for three active neutrinos only, while the R.H.S corresponds to TBM mixing in the 
active part but with extra sterile neutrinos. We remind the reader that sin 2 9 12 = 1/3, sin 2 #13 = 0, 
sin 2 #14 = and sin 2 9 2 4 = on both sides of Eq. (51). Using Eq. (37)-(39) condition (51) leads to 

^ 12-24S 2 , 

s « = Is^ssi (52) 



for a — e, 



2 15 - 72£ 2 3 + 84%, 

bu ~ 14 - 24S 2 3 + 28S& 



3 This is just one example in which we show this confusion. This confusion is expected even in the most general 
case where there might be more sterile neutrinos, with all mixing angles and CP phases non-zero. 
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for a = [l and 

2 27- 96^ + 84^ 

34 27 - 32S 2 2 3 + 28Si 

for a — t. We have plotted Eqs. (52)-(54) in Fig. 4(a), where the solid line corresponds to 
the condition for R e , dashed line for and dot-dashed line for R T . Any point on a given 
line in this figure gives the value of sin 2 #23 and sin 2 #34 such that Eq. (51) is satisfied. For 
instance, the case with three active neutrinos only and with sin 2 2# 2 3 = 0.96 predicts a value of 
R e which can be obtained with TBM mixing in the active sector plus an extra sterile neutrino 
with sin 2 # 34 = 0.245. We note that for the same sin 2 2# 2 3, the values of sin 2 # 34 for the sterile case 
which can also simultaneously satisfy the condition for i? M and R T are different, the corresponding 
values of sin 2 # 34 being 0.21 and 0.11 respectively. This means that in principle if we knew what 
is the true value of sin 2 # 34 , we would be able to distinguish the results for the sterile case with 
that for active neutrinos only with deviation from TBM mixing, if we measured all the three flux 
ratios simultaneously. However, in practice we do not have any precise information on this mixing 
angle. Therefore, this confusion will be hard to solve even if we had measurement on all the three 
flux ratios at neutrino telescopes. 

The second very important aspect evident from Figs. 3 is that the values of the flux ratios 
can easily exceed the range given in Eq. (50). For instance, we see that values of R e > 0.645 is 
possible for sin 2 #23 > 0.4 and sin 2 6> 34 > 0.25. Note that we have kept sin 2 #i 3 = 0, sin 2 9 12 = 1/3 
and all CP phases to 0. In fact for this case, when sin 2 634 = the range of R e corresponds to 
[0.437-0.569], as in Fig. 1 for 5 = 0. We can draw similar conclusions from R^ and R T using the 
middle and right panels of Fig. 3. If the measured flux ratios correspond to such extreme values, 
then this could be a signature for extra sterile neutrinos. 

In Fig. 5 we show the contours of constant flux ratios in the sin 2 # 2 3-sin 2 # 2 4 plane, keeping 
#34 — 0, #14 = 0, #13 = and sin 2 #12 = 1/3. The color scheme is same as in Fig. 3. The R e 
bands increase from 0.436 (lower right-hand corner) to 0.610 (upper left-hand corner) in steps of 
0.017, the i? M bands increase from 0.447 (upper left-hand side) to 0.559 (lower right-hand corner) 
in steps of 0.011, while the R T bands increase from 0.445 (left-hand side) to 0.526 (right-hand 
side) in steps of 0.008. Note that even though we have shown results up to sin 2 6> 24 = 0.2, this 
mixing angle is severely constrained and is expected to be ^ 0.05. We note from the figure that 
R e increases with sin 2 # 2 4, as the denominator in Eq. (43) decreases much faster with sin 2 # 2 4 than 
the numerator. On the other hand decreases with sin 2 # 2 4 because the numerator in Eq. (44) 
decreases with sin 2 # 24 faster than the denominator. For R T the dependence is more complicated. 
The Fig. 4(b) illustrates the confusion that could arise between a three-generation scheme with 
deviation from TBM mixing with the one with sterile neutrinos. In the panel (b) of this figure we 
show the contours corresponding to the condition 

R a (sin 2 # 23 ^ 0.5, sin 2 # 24 = 0) = R a (sin 2 9 23 = 0.5, sin 2 # 24 ^ 0) , (55) 

where a = e, /i or r. As before, the L.H.S of Eq. (55) gives the flux ratios predicted by deviation 
from TBM mixing for three active neutrinos only, while the R.H.S corresponds to TBM mixing 
in the active part but with extra sterile neutrinos. We reiterate that sin 2 #i 2 = 1/3, sin 2 #13 = 0, 
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Figure 5: Contour bands of the flux ratios R e , and R T in the sin 2 #23 — sin 2 #24 plane. The color 
scheme corresponds to flux ratios increasing with red, orange, yellow, green and so on. The R e 
bands increase from 0.436 to 0.610 in steps of 0.017, the R^ bands increase from 0.447 to 0.559 in 
steps of 0.011, while the R T bands increase from 0.445 to 0.526 in steps of 0.008. We have kept 
#14 = = #34 and all other mixing angles at their TBM mixing values. 

sin 2 9 14 = and sin 2 # 34 = on both sides of Eq. (55). From Eqs. (43), (44) and (45) we get the 
conditions 



°24 



4 - 8Sj 3 
9 - 8S 2 3 



(56) 



for a 



q2 
°24 



5 - 24S 2 3 + 285$ 



2:-! 



28S 2 4 3 - 245 2 



(57) 



23 



for a = n and 



S 2 



21 



9 - 325 2 3 + 28S$ 



'21 



1 — 32523 + 28S , 5; 



(58) 



2:; 



for a = r. These conditions are plotted in Fig. 4(b). 

In the general scenario of course we expect both sin 2 # 2 4 and sin 2 # 34 to be non-zero. However, 
we still keep sin 2 # 14 to be zero for simplicity. In Fig. 6 we present three dimensional plots showing 



the flux ratios as a function of sin 2 #24 and sin 2 #34 for sin 2 #23 = 
panel shows i? M and right panel shows R T . We have kept 9 ri 



0. The left panel shows R e , middle 
= and sin 2 #12 = 1/3. 
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Figure 6: Three dimensional plots showing the flux ratios R e (left panel), _R M (middle panel) and 
R T (right panel), as a function of sin 2 #24 and sin 2 #34. The angles 9 13 = = #14, sin 2 9 12 = 1/3 
and sin 2 #23 = 0.5. 

7 Discussions 

In this paper we have considered presence of sterile neutrinos as the only signature of physics 
beyond the standard paradigm. Other forms of new physics could also manifest themselves in 
the observed flux ratios of the ultra high energy neutrinos. In particular, neutrino decay predicts 
spectacular results for the ratios, where one expects the flavor ratio 4:1:1 for the normal 
hierarchy 4 and 0:1:1 for the inverted hierarchy [33]. Authors of [34] considered the effect of 
CPT violation on the flavor ratios, for both stable as well as unstable neutrinos. In [35] prospects 
of probing the existence of pseudo-Dirac neutrinos was discussed. In Table 1 we show the predicted 
flavor ratios for all these new physics scenarios and compare them with that expected for sterile 
neutrinos. We note that the largest deviation for e /0^ comes for the neutrino decay model, 
both with normal (NH) and inverted (IH) hierarchy. The ratio M /0 T however for this model 
stays at 1, just like in the standard neutrino case. This ratio is predicted to be very different 
for the case where we have CPT violation (with and without neutrino decay) and where we have 
sterile neutrinos with large sin 2 # 34 . For sterile neutrinos, this ratio decreases almost linearly with 
sin 2 6*34, the slope being determined by sin 2 6*24. We can see that for very large sin 2 #34, the 4>^/4> T 
ratio predicted by sterile neutrinos is very similar to that expected for antineutrinos with CPT 
violation plus decay and NH. In fact, since M /0 T for sterile neutrinos starts from a large value 
and decreases to a small number. Therefore, for some range of values of sin 2 #34 one would get a 
4>fx/(p T similar to the one predicted by CPT violation with stable neutrinos as well. Therefore, we 
see that there is ample scope of confusing sterile neutrino signatures in ultra high energy flavor 
ratio with that of CPT violation. We also note that the 4> e /4>fi ratio predicted by sterile neutrinos 
could be close to some of the cases for pseudo-Dirac neutrinos. 

So far we have presented all results for ratios defined in terms of the ultra high energy neutrino 
fluxes arriving at Earth. However, one should bear in mind that what is physically relevant as 

4 The authors of [33] get 6:1:1 because they take sin 2 #12 = 1/4. The ratio given above is obtained for TBM 
mixing for which sin 2 612 = 1/3. 
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TV J 11 

Model 


• 2 a 
sin 6*24 


sin 6*34 


rlavor at hartn 


standard 






1:1:1 


Decay (NH) 






/I 1 1 

4:1:1 


Decay (IH) 






0:1:1 


CPT [i($ std + $qpTv)] 






1.07 : 1 : 0.79 


CPT + Decay (NH) [^cptv) only] 






2.43 : 1 : 0.14 


Uri + Decay (lilj [^v^cptv) onl y| 






U.U4 . 1 . o.oU 


Pseudo-Dirac (only 1) 






0.73 : 1 : 1 


Pseudo-Dirac (only 2) 






1:1:1 


Pseudo-Dirac (2 and 3) 






1.43 : 1 : 1 


Sterile 


0.04 


0.0 


1.05 : 1 : 1.04 


Sterile 


0.04 


1.0 


1.05 : 1 : 0.12 


Sterile 


0.1 


0.0 


1.12 : 1 : 1.08 


Sterile 


0.1 


1.0 


1.12 : 1 : 0.32 



Table 1: Comparison of the flavor ratios predicted by the different new physics scenarios considered 
in the literature with the standard three flavor picture (first row) and presence of extra sterile 
neutrinos (bottom rows). Note that for the CPTV+decay models, the flavor ratios correspond 
to antineutrinos only. For pseudo-Dirac neutrinos we show only two cases, one where only the 
oscillations due to the second mass splitting is averaged out and another where oscillations due 
to both second and third splitting are averaged out. We refer the reader to [35] for the complete 
list. Note that we have used mixing angles corresponding for TBM mixing for the ratios given in 
this table. This is why they might differ slightly from the ones given in [33, 34, 35]. 



far as the measurement is concerned, are the ratio of the corresponding number of events in the 
neutrino telescope. The number of events expected in the z/^, v e and v T channels for an assumed 
normalization for the ultra high energy neutrino flux was performed in [32]. While are easiest 
to detect as the resultant muons leave distinct tracks in the detector, the v e are detected by 
observing the electromagnetic showers the electrons create. The efficiency of detecting electron 
events is generally lower than that for muons and it decreases as the neutrino energy increases. 
The threshold for muon events is expected to be about 10 2 GeV while that for electrons would 
be about 10 3 GeV. The v T are detectable through the so-called double-bang and lollipop events 
[36] only at relatively higher energies of about 10 6 GeV. Their detection efficiency would become 
comparable to that for electrons for energies larger than 10 7 GeV [32]. In all our results we had 
tacitly given equal weightage to all three flavors. However from the discussion above, it is clear 
that the behavior of the flavor ratios would change, once we used the event ratios rather than 
the flux ratios, at least for neutrino energies << 10 7 GeV. In fact, for neutrino energies below 
10 6 GeV, we should not include the v T flux at all and work with just the e /0^ ratio. Above 
energies of 10 7 GeV we could in principle work with all three flavors, however the neutrino flux 
falls steeply with energy and we expect only a handful of events at these very high energies. We 
have seen in Table 1 that sterile components cause largest deviation from the standard model 
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prediction for T . To see this we would necessarily need observations at higher energies, where 
statistical uncertainty could be the biggest challenge. A detector much larger than IceCube would 
be probably required for getting an unambiguous signature of sterile species in ultra high energy 
neutrinos. This constraint is applicable to most cases which rely on flavor ratios involving the v T . 
However, we reiterate that even though unambiguous signal for sterile neutrinos in a km 3 detector 
such as IceCube might prove to be difficult, sterile neutrinos will definitely cause confusion with 
the signal for deviation from TBM mixing, and this conclusion is valid for all neutrino telescopes. 

8 Summary and Conclusions 

Neutrino telescopes offer a possible way of constraining the mixing angles. In particular, con- 
siderable interest has been generated of late in using the flux/flavor ratios of ultra high energy 
neutrinos in probing the deviation of the mixing matrix from the TBM mixing ansatz. Since it 
has been observed that the effect of deviating sin 2 9 12 from 1/3 on the flux ratios is not large, we 
kept sin 2 12 fixed at its TBM value of 1/3. Deviation from TBM mixing can then be caused either 
by sin 2 #13 7^ or sin 2 #23 7^ 1/2. We showed that the effect of changing sin 2 #13 from zero is small 
on the flux ratios. We argued that the maximum impact on the flux ratios arise when we change 
the mixing angle sin 2 #23 from 1/2. Therefore, we presented most results by varying the value of 
sin 2 #23 on both sides of its TBM value. 

We presented the expressions for the flux ratios assuming that there was one extra sterile 
neutrino and the mass spectrum followed a 3+1 pattern. Even though we are aware that this 
mass spectrum is now disfavored following the recently declared MiniBooNE results and the mass 
spectrum allowed is the 3+2 scheme, we used the 3+1 scheme as an exemplary case for simplicity. 
Our formalism can be easily extended to the 3+2 scenario, only then one has to contend with many 
more parameters. We showed that even in a very simplified picture where only one active-sterile 
mixing angle sin 2 # 34 (or sin 2 #24) was kept non-zero, it was possible to easily confuse the predicted 
flux ratios for the deviation from TBM case with three active neutrinos with the case where TBM 
holds in the three active sector but we have one extra sterile neutrino. We also showed that for 
very large values of the mixing angle sin 2 #34, which as of now is almost unconstrained, we get 
predictions for the flux ratios that are completely out of the possible range predicted by three 
active neutrinos. If the measured flux ratio really conformed with such extreme values, this would 
then be a smoking gun signal for sterile neutrinos which are heavily mixed with active neutrinos. 
We also presented results where we took sin 2 #24 to be the only non-zero active-sterile mixing angle. 
Finally, we showed results for two non-zero active-sterile mixing angles, sin 2 #34 and sin 2 # 2 4- 

We stress that our results would be valid even if a sterile neutrino is not needed to explain 
the LSND data. Since for ultra high energy neutrinos coming from astrophysical sources, the 
oscillatory terms anyway average out to zero irrespective of the value of Am 2 , a sterile neutrino 
with any value for Am 2 would give such signatures in the neutrino telescopes if the mixing angle 
such as #34, which has not yet been checked in terrestrial experiments, was non-zero and large. 

In conclusion, if the values of the measured flux ratios in the future neutrino telescopes turns 
out be different from their TBM prediction of 1/2, it would still not be a foolproof evidence for 
deviation from TBM mixing. Sterile neutrinos with reasonably small mixing with active neutrinos 
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could mimic a similar response, even when the active sector conforms to TBM mixing. One will 
never be able to tell one scenario from the other by measuring the flux ratios of ultra high energy 
neutrinos, if the measured flux ratios were within the range predicted by the standard picture with 
three active neutrinos. However, if the measured ratios turn out to be clearly outside the range 
predicted by the three active case, we would have a signal for the existence of sterile neutrinos. 
One probably needs a larger detector with sufficient statistics for v T events for an unambiguous 
signal. 
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